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A number of authors have recently pointed out inconsistencies of results obtained with the Huang- 
Yang multipolar pseudo-potential for low-energy scattering [K. Huang and K. C. Yang, Phys. Rev. 
A 105, 767 (1957); later revised in K. Huang, "Statistical Mechanics", (Wiley, New York, 1963)]. 
■ The conceptual validity of their original derivation has been questioned. Here I show that these 

' inconsistencies are rather due to an algebraic mistake made by Huang and Yang. With the corrected 

error, I present the revised version of the multipolar pseudo-potential. 
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5^ ' I. INTRODUCTION 



Most many-body problems require evaluating matrix elements of interparticle interaction in the plane- wave basis. 
For a typical interaction with a short-range repulsive "hard" core such integrals diverge. Nevertheless, this problem 
can be made tractable with so-called pseudo-potential technique, which is usually attributed to Fermi In this 
technique, a suitably chosen pseudo-potential replaces the true interaction; the solution of the Schrodinger equation 
' with the pseudo-potential must reproduce the long-range behavior of the wave-function derived with the original 
(~| . interaction. Since the most fundamental ingredient of the quantum mechanical treatment, the wavefunction, is 
O ' properly recovered, any non-singular property should be well approximated in this approach. 

^ O ^ , In the Huang- Yang (HY) construction [3, 01 j the pseudo-potential Vps is determined as a multipolar expansion over 
delta-function ("lumped") contributions. For a partial wave expansion, ij; (r) = J2i m V'Zm (t) Yim (f), 

> ; Vp^iP (r) = ^ Yi^ if) VI V;™ (r) , (1) 
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where /i is the reduced mass of the interacting pair, k is the conventional collision wavevector, and iji is the phase 
shift for a partial wave I. As to the prefactor /;, the values of the prefactor differ in the original HY paper '3| and in 
"j^ . the Huang's textbook j^, published later. Namely the value of the prefactor from the Huang's textbook is used in the 
literature. Here, by tracing steps in the Huang derivation, I point out that still there remains an algebraic mistake in 
^ the value of the pre- factor fi. The new revised value of the prefactor is 

t3 ■ 
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^ ' where f^^'^'^s jg ^^^^ original (erroneous) pre-factor from Ref.'^. 
• Notice that the s-wave contribution {I — 0) to Vps is not affected by the correction. Since the s waves dominate low- 
' energy collision physics, certain inconsistencies for higher partial waves has not been noticed until very recently 0,0, 
when higher order multipoles became a subject of interest. For example, for identical fermions the s-wave contribution 
vanishes because of the symmetry arguments, and one has to consider p-wave scattering in particular. There are other 
scenarios, e.g., a resonant coupling of c?-waves, when the multipoles beyond I = become relevant. Also the strong 
coupling of higher partial waves to s- waves is a prominent feature for anisotropic (e.g., dipolar) interactions. 

First, Roth and Feldmeier 0| have considered mean- field correction to energies of trapped fermions. They have 
derived their own version of the pseudo-potential aimed at reproducing the energy shifts. By computing the same 
corrections with the HY pseudo-potential, these authors found that each multipolar contribution in the HY pseudo- 
potential must be multiplied by a factor of {21 + !)/(/ + 1) (just as in Eq.©) to bring the computed correction 
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into an agreement with their independent results. These authors concluded that the HY pseudo potential "is not a 
proper effective interaction for a mean field description of dilute quantum gases that goes beyond s-wave interactions" . 
Unfortunately, this statement can be interpreted as if two different versions of the pseudopotential were to be employed: 
one version for continuum and another version for bound-state problems. 

Second, an alternative derivation of the pseudopotential has been presented by Stock et al. @. Instead of delta- 
function lumped at the origin, these authors have proposed to use a shell pseudo-potential, i.e., delta- function placed 
at the surface of spherical shell. In the limit of zero radius of the shell each multipolar term of the HY pseudo-potential 
is obtained. As in Ref. 4J this limit produced multipolar terms differing by a factor of {21 + 1)/{1 + 1) from the HY 
terms. Because of these missing factors, the authors of Ref. Q claimed that there is a "fundamental problem" with 
the Huang- Yang derivation. 

Here I demonstrate that there is no conceptual problem with the Huang- Yang construction. Rather there is an 
algebraic mistake in the derivation Q. I have traced the error to the erroneous application of the Gauss theorem 
in 1^. With the corrected derivation, this additional factor of {21 + 1)/{1 + 1) is fully recovered. There is no need to 
introduce the intermediate (5-shell pseudo-potential as proposed in Ref. In other words, there is a reconciliation 
of the seminal paper by Huang and Yang with the more recent derivations. 

This short paper is organized as follows. In Section m by deriving phase shifts in the Born approximation with 
the HY pseudo-potential I present another demonstration that there is a consistency problem with the original HY 
expressions. In Section IlIII I will point out the algebraic error in Huang- Yang derivation. 



II. INTEGRAL EQUATION FOR THE PHASE SHIFTS AND HUANG- YANG PSEUDO-POTENTIAL 

As discussed in the introduction there is some evidence from the literature 3, Q that there is a difficulty with 
the original HY formula. Below I provide an alternative self-consistency check of the HY pseudopotential based on 
the integral equation for the phase shifts. I arrive at the value of the prefactor which indeed differs from the one 
prescribed by Huang. 

The radial Lippmann-Schwinger equation for a finite-ranged spherically-symmetric potential leads to the following 
implicit equation for the phase shifts 



tan77/ = -J^^/ kr ii{kr)V {r)4)i{r)dr . (3) 



Here (pi (r) is the properly normalized exact solution of the scattering problem. For the HY pseudo-potential this 
solution is by construction 



't'l (^) = a/ TT~r {^^^i i^^) ~ tan??/ kr ni {kr)} , (4) 
h Trk 



with ji{kr) and ni{kr) being the conventional spherical Bessel and Neumann fmictions. At small values of the 
argument, 

. , ^ 1 ; , ^ (2/-1)!! 

Now we substitute the solution Q into the Eq.(|3Jl with the HY pseudo-potential, Eq.(^. The part of the solution 
proportional to the Neumann function vanishes upon differentiation and an intermediate result is 

1 1 
[(2/ -1-1)!!] 

This equation allows us to obtain the value of the prefactor Taking into account (^)^'^^ j,2/+2 _ f^2l + i)! 7-^ one 
arrives at the value of prefactor /; {21 + l)!!/(2'^!) , i.e., differing from the Huang's prefactor by {21 + l)/{l + 1). 
Again, as in Refs. 111,0], we conclude that there is a self-consistency problem with the HY formula for I > multipoles. 
Below, I will trace the steps in Huang's derivation and I will point out the mistake in his derivation. 
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III. TRACING THE ALGEBRAIC ERROR IN HUANG- YANG DERIVATION 

While the original pseudopotential has been introduced in the Huang- Yang paper Q , there are certain mistakes in 
the final equations. These formulae have been revised later by Huang in his textbook Q. Also additional details of 
the derivation are given there. Tracing the steps in the derivation, I found the error in the chain of equations (B.13) 
of Ref. Q . Here Huang integrates over a small spherical volume of infinitesimal radius e 

[ r~ ar ar 



Noticing that V^r' = ^-^^-p^r^ and that for a spherically-symmetric function 

this expression is brought into a form suitable for application of the Green's theorem (the Green's second identity), 

Ie= d^r {r^V^ni{kr) -ni{kr)V^r^} . 



The second Green's identity reads 

Jv Js 



leading to 



while Huang's formula reads 



1^ = j> dS- {r'-Vni {kr) - ni (kr) Vr' 
/Huang ^ ^ rfs • {r^Vni (kr)) . 



The error is here. Huang is missing the second term, — n/ (kr) Vr'. Continuing the chain of equations, we obtain 

le = (f (r''-^ni (kr) - ni (kr) -^r' ] dS = 
Js \ ar dr J 

= (2^-1)!! / (-r^±J- + ^±r^ | dS 
fc'+i Js \ drr'+^ r'+i dr ' 



Namely the second (missed term) provides the additional contribution of I. Finally, 



or Eq.(B.14) of Ref. ^ should read 



4 = 4.(^^(2^ + 1) 



(2^-1)!! A(r) 



Combining this result with the rest of the derivation we arrive at the revised Huang- Yang pseudopotential Eqs. (|1I2|I . 

To summarize, here a revised form of the low-energy multipolar pseudopotential by Huang and Yangj^ ^ has 
been presented. A mistake in the original derivation has been pointed out. The present paper reconciles the seminal 
Huang- Yang construction with more recent results from the literature. 
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